Preconditioning Sparse Nonsymmetric Linear Systems with the Sherman-Morrison Formula
نویسندگان
چکیده
Let Ax = b be a large, sparse, nonsymmetric system of linear equations. A new sparse approximate inverse preconditioning technique for such a class of systems is proposed. We show how the matrix A−1 0 −A−1, where A0 is a nonsingular matrix whose inverse is known or easy to compute, can be factorized in the form UΩV T using the Sherman–Morrison formula. When this factorization process is done incompletely, an approximate factorization may be obtained and used as a preconditioner for Krylov iterative methods. For A0 = sIn, where In is the identity matrix and s is a positive scalar, the existence of the preconditioner for M -matrices is proved. In addition, some numerical experiments obtained for a representative set of matrices are presented. Results show that our approach is comparable with other existing approximate inverse techniques.
منابع مشابه
A Transformation Approach That Makes Spai, Psai and Rsai Procedures Efficient for Large Double Irregular Nonsymmetric Sparse Linear Systems
It has been known that the sparse approximate inverse preconditioning procedures SPAI and PSAI(tol) are costly to construct preconditioners for a large sparse nonsymmetric linear system with the coefficient matrix having at least one relatively dense column. This is also true for SPAI and the recently proposed sparse approximate inverse preconditioning procedure RSAI(tol) procedure when the mat...
متن کاملSome Modifications to Calculate Regression Coefficients in Multiple Linear Regression
In a multiple linear regression model, there are instances where one has to update the regression parameters. In such models as new data become available, by adding one row to the design matrix, the least-squares estimates for the parameters must be updated to reflect the impact of the new data. We will modify two existing methods of calculating regression coefficients in multiple linear regres...
متن کاملDivide and Conquer Low-rank Preconditioning Techniques
This paper presents a preconditioning method based on a recursive multilevel lowrank approximation approach. The basic idea is to recursively divide the problem into two and apply a low-rank approximation to a matrix obtained from the Sherman-Morrison formula. The low-rank approximation may be computed by the partial Singular Value Decomposition (SVD) or it can be approximated by the Lanczos bi...
متن کاملApproximate Inverses and Incomplete Decompositions
In this work we consider the approximate inverse preconditioners AISM based on the Sherman-Morrison formula [1]. We will show some of theoretical properties of the factors of AISM in function of the value of parameter s considered in the decomposition. In particular, we will show that one of the factors hides a scaled factor of the standard incomplete LDU decomposition of the system matrix A. T...
متن کاملLow-rank correction methods for algebraic domain decomposition preconditioners
This paper presents a parallel preconditioning method for distributed sparse linear systems, based on an approximate inverse of the original matrix, that adopts a general framework of distributed sparse matrices and exploits the domain decomposition method and low-rank corrections. The domain decomposition approach decouples the matrix and once inverted, a low-rank approximation is applied by e...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- SIAM J. Scientific Computing
دوره 25 شماره
صفحات -
تاریخ انتشار 2003